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Abstract 

We study the double scaling limit for unitary invariant ensembles of random matrices with 
non analytic potentials and find the asymptotic expansion for the entries of the corresponding 
Jacobi matrix. Our approach is based on the perturbation expansion for the string equations. 
The first order perturbation terms of the Jacobi matrix coefficients are expressed through 
the Hastings-McLeod solution of the Painleve II equation. The limiting reproducing kernel is 
expressed in terms of solutions of the Dirac system of differential equations with a potential 
defined by the first order terms of the expansion. 



1 Introduction 

Unitary invariant ensembles of random matrices or matrix models play a very important role in 
the random matrix theory (RMT) mainly because of its numerous links with another fields of 
mathematics and theoretical physics. An important advantage of these ensembles is that their 
special structure allows to study their limiting eigenvalue distribution with much more details, 
than other models of RMT. 

The matrix model is defined as a set of all n x n Hermitian matrices M with a probability 
distribution 

Pn{M)dM = Z-^ exp{-nTrU(M)}dM, (1.1) 
where is a normalizing constant, U : M ^ M+ is a Holder function satisfying the condition 

U(A)>(2 + e)log(l + |A|). (1.2) 

One of important objects of the investigation in the global regime is the Normalized Counting 
Measure (NCM) of eigenvalues {A^"^}^^;^ of the matrix M. According to [5lll3] the NCM tends 
weakly in probability, as n ^ oo, to the non random limiting measure M known as the Integrated 
Density of States (IDS). The IDS is normalized to unity and it is absolutely continuous if V 
satisfies the Lipshitz condition. The non-negative density p{X) is called the Density of States 
(DOS) of the ensemble. The IDS can be found as a unique solution of a certain variational 
problem [SllTS] which imply, in particularly, that if U'(A) satisfies the Lipshitz conditions on the 
support a of limiting IDS, then p{\) is a solution of the following integral equation 

U'(A)=2/4^, a = suppAA. (1.3) 
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While the global regime depends strongly on the form of V, the local eigenvalue statistics is 
expected to be universal. Denote by Pn{^i, ■■■,^n) the joint eigenvalue probability density. It is 
known (see [15]) that 

n 

Pn{Xi,...Xn) = Z-' n (A,-Afe)2ne-"''^'^\ (1-4) 

l<j<k<n j=l 

where Zn is the respective normalization factor. Let 

pf'\\l,...,\l) = j Pn{Xl, .■.,Xl,Xl+l, ...\n)d\l+l...dXn (1.5) 

be the Ith. marginal distribution density of ()1.4p . Universality of local eigenvalue statistics means 
that if we consider some Aq £ cr, then all marginal distribution densities after a proper rescaling 
(which depends on the behavior of the limiting DOS p(A) near the point A = Aq) tend to some 
universal limits. 

The most known quantity probing universality is the large-n form of the hole probability 

EniAn) = Fn{X\"^ ^ A„, / = 1, n}, (1.6) 

where P„{...} is the probability defined by the distribution (jl.ip . and A„ is an interval of the 
spectral axis, whose order of magnitude is fixed by the condition nAA(A„) ~ 1. For the matrix 
models £^„(A„) can be obtained as the Fredholm determinant of a certain integral operator. 
This structure of the hole probability is a consequence of the structure of marginal densities 
and the latter can be explained by the link of matrix models with orthogonal polynomials 
Pi'^\x), [l = 1, ...) on M associated with the weight e~^^^^\ The link is provided by the formula 

m 

p'^\x,,...,X,) = ^I^-^det\\K^{X,,Xk)\\ik=i, (1-7) 

where 

n 

e:.(a,m) = ^W"^(a)v^;")(m) (1.8) 

1=1 

is known as a reproducing kernel of an orthonormalized system 

V.;")(A) = e-"^W/2p(n)^^^^ ^^^^ ^^_g^ 

m which P/"^(A) is a polynomial of l-th degree with a positive coefficient in front of A'. This 
polynomial is uniquely defined by the orthogonality conditions 

I p/")(A)Pi")(A)e-"^WdA = 6i,m. (1.10) 

Formula (|1.7p gives us 

{n^"^ /coYpi,n{Xo + ti/corC, . . . , Aq + ti/cQV?) = det|/C„(ti,tj)|-j=i, 

where 

ICnitiM) = n^-''cQ^Kn{Xo + ti{con)~\ Aq + t2{con)-''). 

Hence we can reduce the question on the behaviour of the rescaled ith marginal density to the 
question of the existence of the limit of JCn{s, t) for proper chosen 7 and cq. 
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In the bulk case {p{Xo) ^ 0) we choose 7 = 1. Then the hmiting hole probability is the 
Fredholm determinant of the integral operator, defined by the kernel sin7r(ti — t2)/T^{ti — 12) on 
the interval (0, s). This fact for the GUE was established by M. Gaudin in the early 60s [15j . 
The same fact was proved recently in [8] for certain classes of matrix models. In terms of 
the reproducing kernel (|1.8p this result can be formulated as 



lim p-i(Ao)Kn(Ao + ti//o(Ao)n, Aq + t2//o(Ao)n) = ^2) ^ 
"-^00 7r(ti — 12) 

The edge case of local eigenvalue statistics was studied much later even for the GUE [TTl [20] . It 
was found that if we choose 7 = 2/3, then for the edge points Aq = ±a ( a = [—a, o]) the hole 
probability (jl.6p of the GUE in the limit n — > 00 is the Fredholm determinant of the integral 
operator, defined on the interval (0, s) by the kernel 

Ai{ti)Ai'{t2) - Ai'{ti)Ai{t2) 

'^[tl,t2) = ; 



tl-t2 

where Ai{s) is a standard Airy function [Ij. This fact for real analytic potentials in (jl.ip was 
obtained in [9]. In the paper [18] a more simple proof of the edge universality for the same 
class of potentials was given. An important advantage of the method of |18j is that it can be 
generalized to a class of non analytic potentials. 

The case of the critical point universality with 7 = 1/3 was studied first for V{X) = — X'^ 
by using the Riemann-Hilbert approach in [3]. The same method was generalized on a class of 
real analytic potentials in [7] under additional assumptions that the limiting spectrum a consists 
of one interval and the density p{X) behaves like a square root near the edge points and has 
only one critical point inside a (cf. condition C3 below). But the asymptotic behaviour of the 
Jacobi matrix coefficients was not studied. 

In the present paper we find the asymptotic behaviour of the Jacobi matrix coefficients and 
on the basis of this result prove universality near the critical point. We need not to assume that 
V{X) is a real analytic function. Our approach is based on the mathematical version of physical 
ideas proposed in 

Let us state our main conditions. 

CI. The support a of the IDS of the ensemble consists of a single interval: 

a = [-2,2]. 

C2. V{X) is an even real locally Lipshitz function in M and there exists e > such that 
y ( A) e L2 kc] , where = [-2 - e,2 + e]. 

C3. The DOS p{X) has the form 

p(A) = ^A2Po(A)V4-A2, AG [-2,2], (1.11) 

where Po{X) > for X e [-2,2]. 
C4. The function 

n(A) = 2 j log 1^ - X\p{p)dp - V{X) (1.12) 
achieves its maximum if and only if X G [—2,2]. 
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Remark 1. In fact Theorems 1 and 2 below can be proved under the condition V^^\\) € L2[o"e], 
but the proof is more complicated. Since V^^\0) is used in the limiting formulas for TheoremUl 
it is natural to expect that existence of continuous ^(^''(A) in some neighborhood of X = is a 
necessary condition for TheoremUl Thus condition C2 does not look too restrictive. 

Remark 2. It is well known DOS p of the ensemble for a = [—2, 2] has the form 

/o(A) = ^Xa(A)P(A)V4-A2, (1.13) 

where Xo-(A) is the indicator of a and it follows from il.3\) that P{X) can be represented in the 
form 

So condition C3 means that p{X), behaves like square root near the edge points and has the 
second order zero at X = 0. 

Define a semi infinite Jacobi matrix JT^^") , wliose entries ^ = = Ji^^ are defined by 

the recurrent relations 

j;"Vf?i(A) + J,L1#l(A) = AV'J")(A), /GN, Ji")=0, (1.15) 

where V';"^ is defined by (jl.Op . The main result of the paper is 

Theorem 1. Let conditions C1-C4 be fulfilled. Then for any k:\k\< n^/^log^ n 

where q{x) is the Hasting s-McLeod solution of the Painleve II equation 

q"{x) = ,^^xq{x)+2q\x), (1.17) 
which is uniquely defined (see U^ ) by the asymptotic conditions 

hm q{x)=0, lim = , (1.18) 

x^+oo x—>—oo — Xj ' 2i^ (0) 

Po(A) is defined by U.lld . s = sign(l — Jn"^) and remainder terms r^ satisfy the bounds 

\rk\ <c( \k/nf + n-^/A, (1.19) 



where C is some absolute constant. 

To prove universahty of local eigenvalue statistics we study 

/C„(ti, ia) = n-^I^Kn{tin-^/\t2n-^l^). (1.20) 
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Theorem 2. Under conditions C1-C4 for any I G N there exists a weak limit of the marginal 
density il.5]) 

lim (2n2/3)^p(")(2ti/n2/3,...,2t;/n2/3) =det{/C(t,,t,)}J_i, (1.21) 

n-+oo " 

where 

K,[ti,t2) = -^^ , (1.22) 

vr(ri - t2) 

and ^{x,t) = {^Q{x;t),^i{x;t)) is a solution of the Dirac system of equations 

*) = A=[ '_^l)± + [ ) . (1.23) 

with q{x) defined by ^1.17^ - [1.18\) . and ^{x,t) are chosen to satisfy the asymptotic conditions 

lim |*(x;t)|=0, lim |*(x; t)| = 1. (1.24) 

Corollary 1. Under conditions C1-C4 the hole probability il.6\) for A„ = [n^^^'^a,n' 

-1/36] 

converges, as n ^ oo, to the Fredholm determinant of the integral operator defined in [a,b] by 
the kernel M.21\) : 

lim E„([an-l/^ hn-^'^]) = det(I - /C([a, b])). (1.25) 

n— »oo 

The paper is organized as follows. In Section [2] we prove Theorems [T] and [2j The proofs of 
the most of auxiliary results are given in Section [3l Some auxiliary results which have no direct 
links with matrix models (some properties of the Hastings-McLeod solution, bounds for smooth 
functions of Jacobi matrices etc.), are proven in Appendix. 

2 Proofs of Theorems [1], H 

Proof of Theorem [H The main idea of the proof is to use the perturbation expansion of the 
string equations: 

^(n)^,(^H) = -, (2.1) 

n 

(n) 

which we consider as a system of nonlinear equations with respect to the coefficients Jf, . Here 



and below we denote by ^^"^ a semi- infinite Jacobi matrix, defined in (jl.lSp . Relations (|2.1 
can be easily obtained from the identity 

-^Wp("\(A)p("HA))'dA = 



Our first step is the following lemma, proven in Section O 

Lemma 1. Under conditions CI — C3 uniformly in k :\k — n\ < n^/^ 

- 1| < logV4 „^ I + jW^ _ 2| < 

Remark 3. The convergence J^""* — > 1, as n ^ oo and \k — n\ = o(n) without uniform bounds for 
the remainders was proven in [2] under much more weak conditions (V'{\) is a Holder function 
in some neighborhood of the limiting spectrum). 
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Denote j'^'^'^ an infinite Jacobi matrix with constant coefficients 

^(0) _ ^(0) _ -I 

^k,k+l ~ ^k+l,k ~ ^ 

and for any positive N < n define an infinite Jacobi matrix J{N) witli the entries 



Jk 



(n) 



n+k 



1, 



0, 



\k\ < N, 
otherwise. 



(2.2) 



(2.3) 



Proposition!. For any function v{X), whose ith derivative belongs to L2[ue\ (o"£ = [— 2 — e, 2 + 
e]), consider a periodic function v{X) = v{X + 4 + 2e) with the same number of derivatives, and 
such that v{X) = v{X) for \\\ < 2 + e/2. Let also n^/^ >n^m> n^/^ and J{N + M) is defined 
by 12. 3\) . Then uniformly in N , M and \k\ < N for any fixed integer 6 



The proof of the proposition is given in Appendix. 

To estimate the remainder terms of our expansion we define 



(2.4) 



rrik := max 



< max \ \Jj\, 

I lil<|fc|+nl/3/2 I 



Jo- + J. 



i+1 



V2},(|fc|/,)V2 



(2.5) 



Lemma 2. Let v{X) satisfy conditions of Proposition {1\ with i = 5, 6 be any fixed integer and 
\k\ < 3ni/V4. Then 



;(jW; 



n+k,n+k+5 



{J^'^)k,k-,s - c^^Jk + + T!v^ff^JiJ, 



v{J^'^)k,k+5 - c^'^Jk + + + + 



(2.6) 



k ' 



where \r^^^\ < Cm\ and V^^f^''^' and Vi^f^''^^^ satisfy the bounds: 



,(2,fc,<5) 



(3,k,S) 



< C ||x|| ||y| 

< C'll^ll lly|l> 

< C I |x| I I |y| I I \z\\ 



(2.7) 



for any bounded sequences {xk}, {Vk} md {z^}. Here and below \ \x\ \ = max^ \xk\ and means 
the summation over indexes \li\ < \k\ + n^/^ /2 
Moreover, 



1 

2^ 



F('^)(2cos(x/2))e^'^dx, 



with some smooth F^^\X) and for 5 = 1 



4^^ = —/" ?;(2cosx)cosxdx, fW(A) = 2P(A) = - /" 
2vr J^^ vr J_ 



v{X) - v{fi) 



For5 = 



SO) 



0, F(°)(2cos(a;/2)) 



cos^(x/2) 
2^r 



2 (A-^)v^4-/i2 

v{2 cos xi)dxi 
cos^ xi — cos2(x/2) 



(2i 



dfi. (2.9) 



(2.10) 
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The proof of the lemma is given in Section [3l 
Remark 4. If v coinsides with V' for A € then 

IT/ N , r . p(X)dX 

ci = — / i;(2 cos cos x ax = — / ax cosx = 1. (2.11) 

2TrJ_, ^ ^ ttJ^^ y_2 2cos2;-A ^ ' 

Remark 5. // in 112. 9\} P{\) = A^Po(A), then for any Xk 

Y.'^i^-i^i = Y.'^fliixi+i + 2x1 + xz-i), (2.12) 

where 

= - r Po(2cos(x/2))e^'^dx. 

Denote 

Vk = k/{8Po{2)n), (2.13) 
where Po{X) is defined by (jl.lip . Then we represent Jk in the form 

Jfc = (-l)'=Xfc + yfc, or Xk = {-lf{Jk-yk)- (2.14) 

Now, substituting (j2.14|) in (j2.6p and keeping the terms up to the order m| (recah, that by 
definition ([23]) = 0(m|)), we get for (5 = 1 



5,(1) ^ -rpS(-i)'4'^+r^iiW = -sS^^+s5'\ 

2(2) = E'^5f^(-l)'^+'^^h2^^2+2E'^ft'''^(-l)'^^hy/2+0(mT 



(2.15) 



Here and below we denote 

«fe -Xk+i-Xk, % - dfc "fc • l-^-J^oj 

Since by definition (j2.5p < Cml, using bounds (|2.7p . we can write 

= 4E'7'5:'^^(-iy^+'^+2x.d«r<j'^(-i)'^+'H/i-^) 

+2^;^ T: V\lff {{xi, - Xk) - ih - fc)4'^) + Oimt) (2.17) 



+ 2xfc4^)sg + 2x,sg + 0(mt). 
= XkVk T:v\lff{-lY- + 0{ml) = + 0(m|) 



Similarly 



(2.18) 



and 
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Proposition 2. If v{\) = V'{\), as X £ o^ji, then 

= 2(-l)'=Po(0) + 0(n-5/6), = k/n + 0(n-3/2), 

= 1 + Oin'^l^), = 0(n"i/2), (2.19) 

E^^j = (-l)'^ + 0(n-5/6), sl^J = (-l)'=(4Po(0) - 1) + 0(n-i/2). 

//i;(o)(A) = A-iy'(A) /or A G a^/a, i/ien 

Y^vfX = 4Po(2) +0(n-5/6), ^ 2Po(0) + ©(n"^/^). (2.20) 

Substituting (I239D into (12171) - (|2J8]) . we obtain 

V^'(j("))n+fc,n+;c+i = 1 - (-l)'xfc - yfc - j;VS(-l)'4'^ + k/n 

+ xl + 2xfcSg + 2xkyk{-lf + {-lfm{Q) - l)xl + 0{ml). (2.21) 

Using this expression in (j2.ip and keeping the terms up to the order 0{m\), we get 

- E'^i°V-l)^'~'^^f^ + 4Po(0)xE + 8Po(2)xfc2/fc + 2xkT.fl = 0{mt). (2.22) 
Remark 6. If the operator V^^^ has the form i2.12\) with Po{x) > d > 0, then there exists 

= ^ /Vo-i(2cos(x/2))e^('=-')-dx. 
So, if for some Zj we have the system of equations 

Y^V^^liZi = Sk, k€Z, (2.23) 

thcTi 

k.|<El(p(°)),^,||e/|. (2.24) 
I 

We apply this remark to (|2.22p written in the form (j2.23p with Zj. = {—l)'^df'\ where df^ 
are defined by (12.16P with Xk of the form (12.141) for \k\ < n^/^ ^nd Xk = for k > n^^"^. We take 
also 



'4Po(0)x| + 8Po{2)xkyk + 2x,sg + E|,|>iv, ^-/(-l)'^'^ + 0(m|), k < 3nV2/4, 

,E-PS(-l)'t^f^ otherwise. 

(2.25) 

Since P^^\X) has the third derivative (see conditions C2, C3 and representation (jl.l4p . (P^'^^X))^^ 
also does. Hence, using the standard bound for the remainder of the Fourier expansion of the 
functions p(°)(A) and (p(°)(A))~^ (see, e.g. (g^l)), we have for any M 

^ < \{V'^%'\ < CM-'>/\ (2.26) 

|«|>M |«|>M 

Using the first of these bounds and the inequality 

1^(2), 1^(1) _^(i)|<f Cn-V4iogV4^^ \\k\-ny^\>2 
\ Cn"^/^log^/'^n, otherwise 
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which fohows from Lemma [H we get 

'ml + jQgi/2 ^ 



kfcl <C<'n-i/4iogV2^ nV2/2 < <nV2/2 (2.27) 
^Cn~^^^ log^^^ n, otherwise. 

Using this bound in the l.h.s. of (j2.24p and taking into account the second bound in (j2.26p . we 
obtain 

l4'^|<C7m3^[^,/3/4p A;<nV2/3. (2.28) 

Hence, using this bound in (j2.24p and ()2.7p . we get that (see (|2.17p ) is less than 
^(™'t+[ni/3/4])' therefore (I2.22|) can be rewritten as 

- E' T^k-ii-'^ydP + mO)4 + 8Po(2)xfcy, = 0(m^_^[„,/3/4])- (2-29) 
Now subtracting from (|2.29p the same equation written for k := k — 1, we get 

Using Remark [6] by the way described above, we obtain that \d^^^\ < Cm^^^^^-^^.^ for |A;| < 
Hence, writing 

_ ^(2)x-(l) , Y^'^{0), j(2)^ _ ^(2)^(1) , ^,_4 

1/3/2] 



in view of the first relation in (I2.19P , we get from (I2.29P 

4^^ - 22;fc - 2PQ{Q)n ^^ " 1'^*=' - ^^^^^ "^'^ "T'fc+[ni/3/2]- (2-30) 
Lemma 3. There exist C*,L* > such that for any k : > \k\ > L*n^/^ 

mk<C*{\k\/n)^/'^. (2.31) 
Besides, there exist Ci^2,3 si*c/i that for n^^^ < k < k* = [n^^^ log^ n] 

\Xk\ < Cin-i/3e-C^2(fc/nV3)3/2 ^ ^^^4^^_ ^2.32) 

The proof is given in the next section. It is based on the proposition proven in Appendix. 
Proposition 3. Let {xk}\k\<M j satisfy the recursive relations: 

Xk+i-2xk + Xk^i = 2xl + fk, \rk\<£^, l^fcl < ^i- (2.33) 
Then for any \k\ < M - 2Mi with Mi > 2e"V3 

\xk\ < max{e, (2Mf ^ei)^/^}, \xk+i - Xk\ < 4max{e^ (2Mf ^ei)^/^}. (2.34) 
If for \k\ < M 

Xk+i-2xk + Xk-i = fkXk + fk, (2.35) 
with fk> d"^ > 0, then for \k\ < M 

\xk\ <C e-'^l'^-^V,- + |xM|e-'^l*'^-'=l + |x^M|e-'^l^^+'=l. (2.36) 

\3\<M 



Notice that Lemma [2] combined with Lemma [3] give us a useful corollary 
Corollary 2. For any function (j){X) which has two bounded derivatives on [—2 + e, 2 + e] 



n) ' 



n+k,n+k 



(t>{X)dX 



< C(|A;|/n)^/V0(n-i/3) 



Now let us define a continuous function qn{x), which for x G 'Z/n^/'^ coincides with Xk 
and is a linear function for x 'Z/n^^^. Lemma [3] allows us to write (j2.30p as 



(2.37) 



n' 



-2/3 



k 



+ 



k 



+ 0{n 



-1/6n 



(2.38) 



_ni/3 J 2Po(0)ni/3'^" \^ni/3 

where the bound on the remainder is uniform in \k\/n^/'^ < L for any L. We are interested in 
the behaviour of the solution of this discrete equations which satisfies conditions (cf. (I2.3ip and 
(1232)1 ): 

\qn{x)\<C\x\^/^, |g„(x)| < e-^^'^'/2^ as x ^ +oo. (2.39) 

It follows from Lemma [3] that the functions {qnix)}^=i are uniformly bounded and equicontin- 
uous for any bounded interval. Hence, this family is weakly compact in any compact set in M 
and any convergent subsequence converges uniformly to some solution of the Painleve equation 
(jl.l7p . satisfying ()2.39p . Now we need to prove the asymptotic relations (jl.lSp for x — oo. To 
this aim we use Lemma H] below, which describes the behavior of the Stieltjes transform of the 
following densities (cf (fTTT P - ifLS]) ) 

k 
1=1 



Pk,nW = n ^Kk^ri{X,X), gk,n{z) 



Pk,n{X)dX 
z-X ' 



(2.40) 



Lemma 4. For any k : \k\ < n^^^ log^ n gn+k,n{z) can he represented in the form 



9n+k,n{z) 



-{V"{G)z + 



yW(o) 
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-]^X{z) (pi{Q)z* + ^^0(0)^' + Cfc - <5„+fe,„(z) - 4+fc,n(^) 



1/2 



(2.41) 



where X{z) = \/ z"^ — 4 (here and below we choose the branch which behaves like z as z ^ +oo) 
and 

j=0 



Ck = ±n 



-5/3' 



(2.42) 



(ziz corresponds to the sign of k). Moreover, the remainder terms 6n+k,n{z) and 5n+k,n{z) in 
\2.4-l^ for z :\z\ <1 admit the hounds 

i^^+fc,JAi,A2)(Ai-A2)2 



\Sn+k,n{z)\ 



n 



-2 



{x,-zy{X2-zy 



-dXidX2 



< 



C 



1 + 



+ 



(2.43) 



\^n+k,n{^)\ < C 



|2r^+^-2/3 



')+n 



-4/3 



+ 



l og'/' 
nV2 



n \z\ 



v^z 



+ z 
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The proof of the lemma is given in the next section. 

Let us take k = —[Ln^^^] with L big enough. Since it is known (see [12]) that any solution 
of the Painleve II equations which satisfies (12.390 assumes also the bound 

we can conclude that 

0<c,<^+o[{\k\/nf'). (2.45) 

Now let us choose e = n'^/^P^ ""^^^(0) and put in (j2.4ip z = eC,. Then (j2.4ip takes the form 

1 



9n+k,nm = y{0 - ^e2po(o)x(ec) vc^ -e+ck + HO, 

where V is an analytic function, 

< Cfc = P^\0)e-^Ck < — , (2.46) 

(see (1131])), and 

l<^(C)l = Po\0)e~'\^k,nm + knieC) + 0{kn-^)\ < C{1 + |Cp), 
for 1 9(^1 > 1 (see (|2.43|) ). Let b be the smallest root of the quadratic equation 

- LC + Cfc = 0. 

We note, that due to (j2.46p b is real and positive. Consider 

_ PojO) 
2m 

where C consists of two lines 9C = ±1 and 



£ X{eCU{e - h){e -L + b)e-<"/^dC + rL, (2.47) 



' j£ ^(C2-6)(C2-L + 6) - 



Then, using the Cauchy theorem, we get 

I(6,L) = :^^9 j ^J{{exf - 4)(x2 -b){x'^-L + b)e-'''/^dx + fi 
2vr J 

I + I ) V(a;2-6)(x2-L + 6)e-^'/2^x + f l + 0{e) 

J\x\<b J\x\>L-bl 



27r 



= PoiO)Iiib,L)+fL + Oie). (2.49) 

One can prove easily that for large L 

h{b,L)^-CoL^'H^I\ (Co>0). 
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On the other hand, 



I{b,L) = -— I e-" Z^'' limQgn+k,nmdx > 



e^O 



Thus, taking into account (|2.48p 

^1/2^3/2 < ^/|~^| < c"L~^/^ \ck\ = Pl{0)e\L - b)b/A < C{n-^I^L^I^. (2.50) 



The last inequahty combined with (j2.42p . and the bound for the first differences dj^'' imply for 

7 . L + 



-4/30(^1/3^ 



n 



^5/3 Z^^"^ ^1/3' 



2n4/3 



.L-i/6 



n 



-4/3 



((2Po(0)g^(-^)-^ + 0(L-V6)) 



Therefore, 



-L)| = (4Po(0))-^/2l^/2(1 + 0(L-i/2))_ 



But it is easy to show that any bounded for positive x solution of (|1.17p . which possesses the 
above property satisfies also the asymptotic relations 



g„(-L) = s(4Po(0))-i/'L^/2(i + o(L-2)), ~s = signg(O). 



(2.51) 



Hence, we have proved (jl.lSp and now can conclude that qn{x) converge uniformly on any 
compact in M to the Hastings-McLeod solution of (|1.17p . so that 

^n{x) = qn{x) — q{x) ^ 0, as n ^ oo. 

But from (j2.38p we derive that for any x = k/n^^^ and h = n~^^^ we have 

h-^{An {x + h) + An {x-h)-2An (x) ) 
and uniformly in n 



2ql{x) + 2q^{x) + 2g„(x)g(x) + 



2Po(0) 



An{x) + r{x). 



\An{x)\ — > 0, as X ^ iboo. 
Proposition 4. For the Hastings - McLeod solution of there exists 5 > ^ such that 



6g2(x) + 



2Po(0) 



>5\ 



(2.52) 



This proposition allows us to apply the assertion (j2.36p of Proposition [3] to x^ = A{k/n^^^) 
with d = and fk=rk— with from (j2.30p and 

uniformly in k. The assertion of Theorem [1] follows. 
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Proof of Theorem Take some fixed Ci, C2 with 9Ci,2 7^ 0, denote zi^2 = Ci,2'n and 
consider the functions: 

Fn{ClX2) = !{\l-Zi)-\\2-Z2)-\\l-\2fKl{\i,\2)d\id\2, 

J I- (2.53) 

i^i'^(Ci) = n~^'^ {\i-z^r\\2-zi)-\\i-\2fKl{M,\2)d\id\2. 



Changing variables Ai^2 = ^1,2?^ ^^'^1 and using (jl.20p . we get 

Fn(Cl,C2) = j\tl-Cl)~\t2-C2r\h-t2fKl{h,t2)dtldt2. (2.54) 

The proof of Theorem [2] is based on the following proposition: 

Proposition 5. Let the functions Fn and Frp be defined by 112. 53]) and there exists -F(Ci, C2) of 
the form 

i^(Cl,C2) = I l{tl-CirHt2-C2)-Htl-t2?Htut2)dtidt2 

with $(ti,t2) bounded uniformly in each compact in and such that uniformly in 9Ci,2 ^ 1 

|Fn(Ci, C2) - F(Ci, C2)| < Cil + |CP)n"i/^ (2.55) 
Let also uniformly in Im( > Sn = (logn)~^/^ and '^C, varying in any compact L C M 

|F«(C)|<C(L). (2.56) 

Then for any intervals Ii,l2 C M 

lim / dti / dt2lCl{ti,t2) = / dti / dt2<^{ti,t2). 
Proof of Proposition 0. Notice that 

where Sk,n{z) is defined in (|2.43p . Therefore, using the bound (12.56P in (j2.41l) . we get for any 

l5n,n(n-i/=^C) + 9nA-n~'^'0\ < Cn-'^'Hd^ + 1), (2.57) 

where C does not depend on n and C- On the other hand, taking z = n^^^^{a + ie), we have for 
any e > e„ 



/C„(t,t)dt < 2e' [ = 2en2/3c>5„,„(.) 

-a\<e J (t-a)^ + e^ 

< Snn^^H'^9n,ni^) + %n,n(-^)) < eCia^ + 1). (2.58) 

Take the integral 



/ dCi [ dC2{Fn{Ci,C2)-F{Ci,C2))e- 



(Ci-ai)V2aig-{C2-a2)V2<T2^ 



Using the Cauchy theorem, we get that for any ai^2 > 0, 01^2 G ^ 

{ti-t2f {Kl{ti,t2) - $(tl,t2)) e-(*i-'^i)'/2<xig-(fe-a2)V2a2^^^^^^ 
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with C, depending on ai, 02, cJi, a2, but independent of n. This imphes that for any Lipshitz fi 
and /2 with a compact support 

Urn [{h-t2f {lCl{h,t2)-Hh,t2)) h{ti)f2{t2)dtidt2 = 0. (2.59) 

n— >oo J ' 

For any small enough e denote by /^^"^^ a Lipshitz function which coincides with the indicator 
Xii of Ii = (ai, 61) inside this interval, equals to zero outside of (ai — e,bi + e) and is linear in 
(oi — e, ai), (61, 61 + e). Let f[ be a similar function for the interval (ai + e,bi — e) and 72^*^^ 
be similar functions for I2. Denote also 

</'ei(il,*2) = {h - t2)~^l|ti-t2|>ei + ^r^l|ii-t2|<ei • 

Then, evidently 

CpeAh,t2)fi-'\h)ft\t2) < Mtut2)xlAh)Xh{t2) < (t^sAhM) f^^'\tl) i^'\t2) ■ 

Integrate this inequality with {ti — t2)^^n(^i) ^2)) and take the limits n — > 00 and then e — > 0. 
We obtain 

/ dtidt2^{ti,t2) - 0{ei) < Yim I dtidt2lCl{ti,t2) 

+ lim / dtidt2lCl{ti,t2){e^'^ih-t2f -l)lit,-t,\<e^< [ dtidt2'^ih,t2)+0{ei). 

But using the inequality /C^(ti,i2) < ^n(ii) ii)^n(i2, ^2) and then (I2.58p . we obtain that the 
second limit is 0(ei). Then, taking the limit ei — > we get the assertion of Proposition [5l 

Let us check that in our case conditions ()2.55p and (|2.56p are satisfied. Using the Christoffel- 
Darboux formula, it is easy to derive from (I2.53P that 

..... . ,(n).2 / (^i"^(Al)V^a(A2)-^i"^(A2)ri"jl(Al))^ , 

"^^^^^'^^^ = ("^^ ) y {X, - z,){X2 - Z2) ^^^^^ 

{Z2) + Rn,n{z2)Rn+l,n+l{zi) (2-60) 
- 2Rn,n+lizi)Rn-l,niz2)], 

Pn\Cl) = 2{J^^f[Rn^n{zi)Rn+l,n+l{^l) " Rn,n+l{^l)Rn,n+li^l)]^ 



where 



Rk,m{z) = I ^ _ ^ dX (2.61) 



is the resolvent of {R = {z - jW)"!). 

Proposition 6. Xei J be an arbitrary Jacobi matrix , with \J'j,j+i\ < A, for all j such that 
\j -k\<M. Consider lZ{z) = {z - J)'^ with < Ai. Then 

\n,,{z)\ < gle-^^l^^ll'^-^-l + ^e-^^l^^l^^ (2.62) 

where Cj, C2 > depend only on A and Ai. 
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Let us choose = [n^/^ log^ n]. By ()2.62p . for a = 0, 1 and 9z > n~^/^e„ 

fc:|fc-n|>Ar 

Therefore for a, /? = 0, 1 

K+a,n+/3(^) = E ^n+a,fc(^)i?n+Afc(^) +0(e-^l°^'''"). (2.63) 
|fc|<Af 

Consider the matrix J'^'^<'^^) whose entries coincide with that of ^T^") with the only exception 
Jn±2N]n±2N+i = 0" Consider 

We need also a simple observation, following from (j2.62p and the resolvent identity 

_ i?{0) = i?(0)(M(i) - M(o))i?(i), i?(i-0) = _ m(1'0))-i. (2.64) 

Remark 7. // in ^/le resolvent identity we take M(°) = J'^") and M^^) = j("'2^) iV^ = 
+ v}^'^ log^ n, t/ien ^2. 62\) gives us for any z : \Qz\ > n^^l'^ log^"*^/^ n for any fixed 5 

\Rn+k,n+,^s{z) - R^^X+k+s\ < C^e-^^l°s'''^ (2.65) 

Let us study first the case when in (jl.l6p s = 1. Consider the Dirac operator A defined in 
L2{R) X L2{R) by the differential expression (fT^Sil - ifLTSl) . Let TZa,fiix,y;C) {a,(3 = 0,1) be the 
kernel of the operator 7^(C) = (C-2^)-^ It means that the coefficients TZa^p{x, y\ C) satisfy the 
equations 

2-^^1,0(2;, y; C) + 2g(x)7?.i.o(a;,y;C) - C^o,o(2;,y;C) = S{x-y) 

d 

-2-l-■7^o,l(a;,y;C) + 2g(x)7^o,l(a;,y;C) - C7^i,i(^'y;C) = K^-v) 

(2.66) 

-2—TZofl{x,y-Q+2q{x)'RQ,Q{x,y;C)-C,'Rifi{x,y;C) = 
2^7^l,l(x,y;C) + 2g(x)7^l,l(2;,y;C)-C7^o,l(2;>y;C) = 

Here 5{x) is the Dirac (5-function and, e.g., the first equation means that the l.h.s. is equal to 
zero, as x 7^ y and 27^.1^0(2:^ + 0, x) — 2TZifl{x — 0, x) = 1. 
Consider the (4A^ + 1) x (4A^ + 1) matrix with entries 

-( ^\ik+^)v /^ 2A; + a 2m + /? \ 

^n+2k+a,n+2m+f5 — \ ^) '•*-Q!,/3 I ^^^3 , ^-yj^ ) ' V'^-^'l 

where — < k^m < N — \, a, (3 = Let us prove, that for \k\, \m\ < N 

i?S,n+^(Cn~^/^) - K+.,n+™| < Cn-''\l + \C?). (2.68) 
Using equations (|2.66p . one can check that for all —2N < k,m < 2N, a, /3 = 0, 1 

{{Z - J("'2^))i?*)„+fc,„+^ = 5k,m + 4,m + a^rk^5k,2N + agHk,-2N (2.69) 
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where the remainder terms could be esthnated as follows: 

|4,„^| < Cn-2/3 ^ (l + |Cp + g2(x))|7^,,^(x,y;C)| 

a,/3=0,l 

\dk,k\ < Cn-^l^ {^ + \C\+q{x))\naAx,y;C)\ 



_fc + a_ m + e(a,/3,fc,m) 

nl/3'«- „l/3 



(2.70) 



fc + Q _ fc + fl(a,^.fc,fc) 

■riTa'S' 1173 



with /3 fc ml ^2. The components of the vectors a^^^ and a^-^^ satisfy the bounds 



I^') I 



„| < |7^o,a((2iV + l)n^l/^(2m + a)n-l/3;C)l +Cn-i/3(|^| + |log2 



n 



< |7^o,a((-2iV - l)n-i/3, (2m + a)n-^/^;C)\ + C7n-^/3(|^| ^ 1 1^^2 ^j^^ 



(2.71) 



n+i,n+j — dn+i,n+j, Dn+i,n+j — O-j ^ Si,2N + O-j ' 5i-2N ■ 



Let us define the (4iV + 1) x (4A^ + 1) matrices 
D 

Then ()2.69p can be rewritten as 

_ j(n,27V))^* =/ + i? + i) 



,(2); 



(2.72) 



Using the trivial bound for the norm of any matrix \\A\\^ < TrAA* and the bound for the 
resolvent of the Dirac operator (see [1] 



\na,p{x,y;C)fdy < C\QC\-\\q{x)\ + \C\), 
we obtain from the first and the second line of (j2.70p that for > e„ 

ll^^ll < C{\C\ + l)n~^/3log=^n < C7(|C| + l)n-^/^. 
On the other hand, it follows from (|2.72p that for any —2N <k,m< 2N 



n+k,n+m 

and so we can write 



n+m 



|4T - {R*{i + D)-XA < \\D\\ . ll^ll^q^p^^^^^^ 

But for any x G M^^+i, 

(/ + I))x = X + e„+2Af(x, a^^)) + e„„2Af (x, a^^^), 
where {^k}k=n-2N ^ standard basis in IR^^+^. So we get 



\\{I + D)-^\\<\\A-^\\ + l, A 



1 + a(^) a(^) 

^ ^ "'n+2N-l "'n~2N 



^n+2N-l 



^n~2N 



(2.73) 



(2.74) 



(2.75) 



(2.76) 



(2.77) 



The off diagonal entries of A tend to zero, because they are expressed in terms of the resolvent 
T^a,i3{x,y; C) with |x — y| ~ log^n (see (j2.7ip ). a'^^l2N-i ~^ 0' because q{x) ~ C\/—x, as 

1/4, because if the potential tends to zero fast enough, then the 



-DO. And la, 



(2) I 
'n+2N\ 
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resolvent near the diagonal coincides asymptotically with that for q = (see [14J). Hence, using 
(|275]) and (^7!7\\ . we get 

R^,T -{R*{I + D)-%m\<Cn-'^\l + \C\'). (2.78) 
But from (I2.76|) we derive 

Hence, we obtain from (I2.78|) and ()2.62p 

l4T -^^nl < Cn-y' + Wil + Dr'\\i\K^,j,\ + < Cn-'/' (2.79) 

and we have proved (j2.55p with 

F(Ci, C2) = (7^o,o(0, 0; Ci)^i,i(0, 0; C2) + 7^o,o(0, 0; C2)7^l,l(0, 0, Ci) 

- 27^o,l(0, + 0; Cl)7^l,o(0 + 0, 0; C2)), (2.80) 

where we denote 7^0,1(0, 0+0, Ci) = linix^+o ^o,i(0, a^, Ci)- But according to the spectral theorem 
(see [H]), 

TZa,p{x,y;C) = J ^3^^ at, (2.81) 

where *(x;t) = (^'o(0; t), ^i(0; t)) is the solution of the Dirac system (jl.23p . satisfying asymp- 
totic conditions (jl.24p . The last two relations and the formula of the inverse Stieltjes transform 
yield 

^.(ti, t2) = (2^)-' (^-1(0; ti/2)^o(0; t2/2) - ^-0(0; ^2/2)^-1(0; . (2.82) 

Moreover, similarly to (|2.79p we obtain 



|A;-n|<Ar 



= -^no,oiO,0;C)+0{\\D\\ . I^Cr')- (2.83) 



Using the representation (j2.8ip and taking into account that ^oiix',t) are smooth function with 
respect to t, according to the standard theory of the Cauchy type integrals (see [16]) we get 
that the derivative in the r.h.s. of (j2.83p is uniformly bounded up to the real line. Therefore 
we obtain (j2.56p and then, on the basis of Proposition [5]) , obtain the assertion of Theorem [2] for 
1 = 2. For the other I we study by the same way 

FniCl,---Xl) = [HiU - drHh - t2)...{tl - h)ICniti,t2)...ICniti,h)dti...dti 

i=l 

Now, notice that the {^Q{x,t),^i{x,t)) (-^'i(x, t), ^'o(x, t)) gives us the solution of (|1.23p 
with potential qi{x) = —q{x) but does not change the expression (jl.22p . This completes the 
proof of Theorem [2] 

To prove Corollary 1 we split the expansion for the Fredholm determinant in two parts: with 
m < N and m > N (m is the number of variables in the correspondent determinant). Using 
the Hadamard bound for determinants with m > N and then (j2.58p it is easy to see that the 
second sum possesses the bound /N\. Hence using Theorem [2] we can take the limit n — > 00 
in the first sum and then take the limit N 00. Relation (jl.25p follows. 
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3 Auxiliary results 

Proof of Lemma [H We introduce an eigenvalue distribution which is more general than (jl.4p , 
making different the number of variable and the large parameter in front of V in the exponent 
of the r.h.s of (fO|l : 

k 

Pk,n{^i,...Xk) = Z^l n (Ai-A™)2expne-"^^'^\ (3.1) 

l<j<m<k j=l 

where „ is the normalizing factor. For k = n this probability distribution density coincides 
with (O). Let 



Pk,ni^l) = J dX2...dXkPk,niM, ■■■h), ^fc,n(Al, A2) = y dX^...d\kPk,n{^l, ■■■Xk) (3.2) 

be the first and the second marginal densities of (jS.ip . By the standard argument [l5j we obtain 

Pk,n{X) = k-^Kk,n{X,X), 

Pk,n{X,^l) = yj—-^[Kk^n{X,X)Kk,n{p,p)-Kln{X,^l)], 

where Kk^n{X,p) is defined in ()2.40p . Remark also that 

Tl 

Pk,n{X) = -j^Pk,n{X), 

where pk^n is defined in (I2.40p . Taking any twice differentiable and vanishing outside cJ2e function 
0(A) and integrating by parts with respect to V , we come to the identity 

V'{X)h,nW^WdX = - [ pk,nW<t>'{X)dX + 2^ [ pk^^{X,^)p^dXdfi. (3.4) 

n J n J X- fi 

The symmetry property Pk,niX, p) = Pk,n{p, A) of (j3.2p implies 

Pk,n{^,p)-p—^dXdfi = - I pk,n{X,p) f^^^ dXdfj.. 
X- p J X- fi 

This allows us to rewrite (13.41) in the form 



V'{X)pk,n{X)HX)dX = - [ pk,^{X)^'{\)dX + ^ I /5fc,„(A,M)^%-^dAd/.. 

n J n J A — p, 



Now, using (|3.3p and the fact that 



j Kln{X,p)dp = Kk,n{X,X), 



we can rewrite the last equation as 

m - <p{p) 



X — p 
where we denote 



Pk,n{X)pk,n{p)dXdp - j V'{X)pk,niX)(t>{X)dX + 6k,n{(p) = 0, (3.5) 



6k,ni'^) = ^ J [<P\X)+<t>\p)-2^^-^^Kl,{X,p)dXdp. 



(3.6) 
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Subtracting from ()3.5p the relation obtained from ()3.5I) by the replacement k ^ {k — 1) and 
multiplying the difference by n, we obtain: 



2 / ^^^^-|^p(^)[V^l")(A)]2dAdM-| y'(A)0(A)[^l")(A)]2dA 



+ 42('^) + S('^)=0, (3.7) 



where 



= / - p(^))[Vi"^(A)]^dAd/. - 1 1 <^'(A)[^l")(A)]2dA. 

By Schwartz inequality 

■ ( / (V'l"^(A))2(4")(^))2dAd^y^' < ^||<^"||o (3.9) 

iC(^)i < + ^ii^'iio < cfii^'iif + ii^'iio^ 



where the symbols ||...||o and ||...||2 denotes the supremum and the L2-norm on a^. Here we 
have used the result of [5], valid for any smooth function (/>(/u) defined on 



<C||,/.'||f ||</.||f n-V2iogV2„^ (3.10) 

where the symbol ||...||2 denotes the L2-norm on a^. 

Now we are going to use (jl.Sp in the second integral in the r.h.s. of ()3.7p . But since this 
representation is valid only for A E [—2, 2] we need to restrict the integrals in (j3.7p by some 

= [—2 — e,2 + e] with some small e > 0. To this aim we use 

Proposition 7. Consider any unitary invariant ensemble of the form U.l\) and assume that 
V{\) possess two bounded derivatives in some neighborhood of the support a of the density 
of states p. Let also a consist of a finite number of intervals, p{X) satisfy condition C4 and 
p{X) ~ C{a*)\X — a*\^^'^ near any edge point a* of a. 

Then there exist absolute constants C,Co,eo > such that for any positive Con~^^^logn < 
£ < Eq and for any integer A; : < n + n^^'^ the bounds hold: 

nCe f ^,/,(")/\\\2j\ ^ ^-nCe 



PkAX)dX<e~^^', / {^r{X)Yd\<e-^^\ (3.11) 

R\(T£ JR\o-e 

This proposition was proved in [3]. It allows us to restrict the integration in the first three 
integrals of (j3.7p by with e = CQn~^/'^ logn. Now we can use (ll.3p . The error, which appear 
because of this replacement is of the order 0{e), because V'{X) is a smooth function in a^- 
Hence, (13.71) can be rewritten in the form 



2 / (V'i"^(A))^dA ii^p(^)d^ = 4«)('A) + C^(<^) + 0(ll'^l|on-^/2logn). (3.12) 
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Take (/)(A) = Pq^{X){z - X)-^ and substitute in (IXT21) . Then, according to (jl.lSp . we get 



2 (V'r(A))2dA 



df^ = sf^iz) + ~6f^{z) + Oi\Qz\~^n-^/Hogn), (3.13) 



-2 

where 5^^^ (z) and 6^/^^ (z) have the form (j3.8p and due to (|3.9p satisfy the bound 



IC(-)I < 



I 19' 



Thus, using the fact that 



Try (z-^)(A-^) 



dfj. 



\6^^Hz)\<^ + ^^^ 



— ^ ^ {z^ + zX + X^) + 2, 

z — A 



(3.14) 



we get from p.l3p 

RkM^) = {z' + ak + 45 (^) + S(^) + 0(|9z|-in-i/2 logn) 
where Rk^k{z) is defined in (j2.6ip and we denote 

«fc= /" A2[V'(")(A)]2dA-2. 



1 



2:2^/^234' 



(3.15) 



(3.16) 



Let us assume that > Cn log^^^ n with C big enough. Then, using the bound (j3.14p and 

1 1/2 

the Rouchet theorem, we get that Rk,k{z) has a root in the circle of radius -^o-jJ centered in the 
point ia^ . But, by definition (j2.6ip . 



9iife,fc(z)9z < 0, 



(3.17) 



so Rk,k{z) cannot have zeros, when Qz ^ and therefore we get < Cn ^/^log^^^ n. Similarly, 
if we assume that Ofc < — Cn~^/'^ log"'^'^^ n we get that QRk,k{^\o.k\^^'^e'^'^^^) > 0, which also 
contradict to ()3.17p . Thus, we obtain that 



\ak\ < CrT^I^Xo^l'^n 

From p.l6p and (|3.15p we find 



(3.18) 



A2(^(")(A))2dA = 2 + afc 



A4(V'f (A))2dA = — ^ CRkAQdC, 



(3.19) 



6 + 2afc + 0(n-i/2 logn). 



Using here the first equation for := zt 1 to express (J^!j_\)^ and through {J)!^'Y , we 
obtain 

I o I 1 1/2 

V 0{n ' logn) . (3.20) 



7-(") \2 



7(")\2 



- 1 + Y + 4 ± 



Combining this relation with (j3.18p , we get the first statement of Lemma [H The second state- 
ment follows from the first one and the first equation of ()3.19p . 
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□ 

Proof of Lemma 

Choose M = cn^^^, where the constant c is small enough to provide the condition 

dC2C < Ci/7, (3.21) 



where Ci and C2 are the constants from (j4.5p and d = 7r(2 + e) ^. This condition and (j4.5p 
guarantee that for any 1,1' : |/ — l'\ > and any j : \ j\ < M, \t\ < 1 

Applying (14. 3p three times we get p.6p with 



k,k+S 



-1+S2 = 1 ^ 

» A/ 

Pgif ) = / dsids2dss ^ji^^df (e^i'^^i^^°'i?(^i)e^i'^^2^^°'E('2)e 



» M 

',k+5 



k,k+5 ' 



E«.=i ,=_M (3.23) 

\ J k,i 



'k 



M 



^ / dsi...ds^Y Vj{ijd)^ 

«i,...,U-^Es»=l j=-M 

\,ijdsij(0) J^ijdS2j'-0) J^ijdssJ'-O) J^ijdS4j<-°^ Jf,ijdS5{j(0)^J) 



+ / dsids2 ^ Vjiijd) U^'^''^^"^ J{e 

Jsi+S2 = l 

where we denote by E^''^ a matrix with entries: 



k,k+S 



ijdsij(O) j(^ijds2j^°'> _ ^ijds2{j(°'>+J)- 



k,k+S ' 



-^fc,m ~ ^k,lSm,l+l + Sk,l+l^m,l- 

Using the Schwartz inequality, we have 

1/2 / \ 1/2 



Hence, using again the Schwartz inequality, we obtain 



< m|d^ b'l^l^il + "^kd bll^il ^ + C"^fc^"^^^ < C'"^t., (3.24) 

\j\<M \j\>M 



where the last inequality is valid because of the choice of M and (|2.5p . 
To obtain ()2.7p we use the representation (see jl]): 



l^^ijdsjW-^^ ^ = J_r ^ijdscosx^i{k-l)x^^ ^ Jk-iUds), (3.25) 
27r J_^ 
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where Jfc(s) is the Bessel function. But it is well known (see, e.g. [Ij) that the Bessel functions 
satisfy the following recurrent relations: 

k^kis) = ^ i Jfe+i(s) + Jfc-i(s) ) • 



2 

Thus, e.g., the first sum in (j2.7p can be expressed via the terms 



Vj{ijdf I dsids2dsz{djsi + ai){djsi + a^) 

\3\<M •^E«.=l 

• ^ Jfe_i,+„3(djSi)i/^J/^_/2+Q4(c|7S2)y/2J/2-A;+a5(2jS3), 

where qi, . . . 05 can take the values 0, ±1, ±(5 + 1). It is easy to see that any of these sums can 
be written in the form: 

where evidently 

< max|xfc|, < max|yfc|. 

Hence, similarly to (j3.24p we obtain 

T!'Plf:t''\h-kfxi,m2 <C\\x\\\\y\\ \j\%\<C\\x\\\\y\\. 

\j\<M 

The other inequalities in (j2.7p can be proved similarly. 

We are left to prove (j2.9p . Due to representations (|3.23p and ()3.25p . we derive that V^^"^ can 
be represented in the form (j2.8p with 



Using pin and ([3:251) we get 

Jo f J — Tl J — TT 

• exp{2zjd[si cos xi + (1 — si) cos X2Wdxidx2 
(1) 1 /■'^ v(2 cos xi) — i;(2 cos(xi — x)) 



27r 



/'^ v(2 cos xi ) — i;(2 cos(xi — x) ) , , , 
^ — — -il + cos 2x1 - X dxi 
cosxi — cos(xi — x) 

r 1 + 

/ i;(2cosxi) 

P(2cos(x/2)) +P(-2cos(x/2)) 



(1) 1 N 1 + cos(2xi - x) 1 + cos(2xi + x) , 

1 H / v(2 cos xi) ^ H dxi 

zn \cos xi — cos(xi — xj cos xi — cos(xi + x) / 



(3.26) 

Representation (j2.10p can be obtained similarly. Lemma [2] is proven. 

Proof of Proposition\^ Let us remark first that all limiting expression in the r.h.s. of (|2.19l) 
and (j2.20p correspond to infinite sums over j in the definitions (j3.23p and infinite sums with 
respect to all /j. The estimates for the remainder terms, which appears because of the restriction 
of summation in (j3.23p over \j\ < M, were obtained already in the proof of Lemma [2j And the 
remainders, which appear because of the replacement of infinite sums by sums over < N^, 
can be estimated by 0{e~^^'^^ ^/12^ ^j^g I^q (j3.22p . Thus we are left to compute infinite over li 
sums for 'p^'^f'^'^ and p^^f'^'^ 
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The first relation in ()2.19p follows immediately from ()2.8p and ()3.26p . To obtain the others 
let us consider an infinite Jacobi matrix J^'"^ with J^'^+i ~ "^fc+i k ~ define 



27r(a2 - 62) 

-ifh 



y^(A)signA „ L,^ d\ 



+ 



27r(a2 - 62) 



(4a2 - A2)i/2 
y-(A)signA |^,_^^,j^^, dA. (3.27) 



It is easy to see, e.g., that 



2 362 



1 r Ay'(A)dA 1 r V'{X)dX 



a=l,b=0 



2vr \/4- A2 TT y_2 AV4 - A2 



1. 



Here we have used (]2.1ip for the first integral and (I1.14[) the second. Similarly 



2dadb 



Vkia,b) 



(-!)'■ 



a=l,b=0 



To compute the sum for 'Pi'^f^']]^ let us observe that 



6(963 



6=0 



(-1)^ 1 

2 962 (1 _ 52) 



7(6) 



6=0 



where 



m = ^ j{V'{X) - Ay"(0))signA^^i-^dA 



^ V'iX) -^SnX^^-^dX - V"mi - 6^). 



Differentiating this expression, one can easily get the expression of (|2.19p . 

To prove the last relation in (j2.19p we use the symmetry arguments. Indeed, according to 

h{kM) = E(^ft'/^ = Y.''jmf I dsrds^ds, 

1 -■ Si +So + S.l = l 



where 



Usi {k - h)fs2,s3 {k-h)+ {k-h - l)/s2,s3 (^1 -ki-1) 



usdk -l) = (e*^'^^^^'°')fc_/, /.„.3(/ -k) = (-l)'(e*^'^^^^*"' jWe*^'^^«^'°^,fc. 



Since both Us^ {k — I) and fs2,s3 (l — k) are even functions with respect to {l — k), after integration 
with respect to si, S2, S3 we get that 

h{k, h) = h{k - h) = h{h -k)^Y, ^(^' - ^1) = 0- 

h 
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To prove in:2Qh we define similarly to IK27\\ 

1 f iV'iX)-XV"m signX 

Then 

(2,M), iy,+/. _ 1^^0.1 _ 2 /-^ (r(A) - XVmdX _ 

^ ~ 2db^^'^^^'^\^,- 7rJ_, A3X(A) -2^0(0). 

Proof of Lemma\^ Relation (j2.30p can be written as 

df = 2xl + fk, h = Xkk{2Po{0)n)-^ +rk, \h\ < C4mk\k\/n + mt), (3.28) 

where C=k is independent of N, n and we always can choose C* > 1. If m^. < for all k > n^^^, 
then (j2.3ip is fulfilled. If m^, > k~^ for some A; > Tn^/^, we can apply Proposition [3] to {xj}y^<^M} 
with M = k, Ml = [n^/3/2], = (?tia;+2Mi (^ + 2Mi)/n + ?Ti-fc_,_2A/i ) ' = i^k+2Mi, because 
Ml > 2/3e"^ Then, since 

2eiMf ^ = 8mfc+2A/in"2/^ < C^mk+2MAk + 2Mi)/n < e^ 

we obtain by (j2.34p that 8e^ = 8C* {rhk+2Mx{k + 2Mi)/n + ?ti|_^2Mi) — ""^fc- Therefore at least 
one of the following inequalities holds 

8amfc+2Afi(fc + 2Mi)/n > ml/2 V 8C*m|+2A/i > (3.29) 

Since according to Lemma [U | mfc_|_2A^i I < Cn^ log^''^ n the second inequality yields 

rhk+2Mi > 2Thk (3.30) 

If the second inequality in (j3.29p is false, then the first one holds. Write it as 

rhk+2M, > {lQC,)-^mk [rhln/k] \k/{k + n^/^)] . (3.31) 
Assume that for some k > 7n^^^ 

mln/k > 32C^. (3.32) 

Then dSDimphes fjjmi and 

^l+2M,n/{k + 2Mi) > 4 [mln/k] ■ {k/{k + 2Mi)) > 32C,. 
Hence, we can repeat this procedure / times with / = [logn]. Then we obtain the inequality 

rn'k+[\ogn]A4i > 2[^°§"lmfe, 

which contradicts to Lemma [TJ Thus, (j3.32p is false and we have proved (j2.3ip . 

To prove (j2.32p take any ko > n^/^ denote Xk = Xk^2ko smd, taking into account (|2.3Up . 
apply (|236]) with M = ko. Then since fk > {ko/2Po{0)n) we obtain ([2:32]) . 

Proof of Lemma 
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Substituting in ()3.5p (/)(A) = {z — X) ^ we get easily the equation 

/V'(X) 
YrxPn+k,nWd\ = 5n+k,n{z) (3.33) 



with 5n+k,n{z) of the form (cf. (13. 6p ) 

(5„+fc,„(z) = ""^^ ((ii2)„+fc,„+fc(^;)(^^)n+fc-l,n+fc-l(^;) - (i?^)^+fc,„+fc_l(^;)) , (3.34) 



where (i?^)^^^. = {z - J^"))^]. Here we have used the Christoffel-Darboux formula in the 
numerator of the first integral in (j2.43p . Let us transform 

^A) , , 3^^'H0)\ fV'{\) 



z-X 



Pn+k,nWdX = gn+k,n[z) I zV (0) + Z I - / ^ Pn+fc,n(A)dA 



[V'{X)-xv"{Q) /■ y-(A)-Ay--(o)-iA3yW(o) 

-Pn+A:,n(,AjaA - Z / -j- — Pn+k,n\X)aX 



A3 ™,-V'V-- y A4(z-A) 



= 9n^kA^) \ zV'iO) + z^^] - - z^cgi - /ca(.). (3.35) 

Denote 



Taking the limit n — > oo in (j3.33p and using (jl.lSp . we get for any A G [—2, 2] 

X'P^{X){X' - 4) = [V'{X)f - 4Q(A) Q{X) = \ {[V'{X)f + X'Pi{X){A - X^)) . (3.36) 
Therefore, denoting v^^\X) = y'(A)A"\ we get 



~{0) 
Cn 



/ t;(°)(A)(p„,„(A) -p(A))dA, cu = in"^ ^ t;(°)(j("))„±,, 



(3.37) 



Here and below in the proof of Lemma H] the sign it corresponds to the sign of k. Repeating the 
argument of Lemma [2] for the function v^'^^X), we obtain 

\h\,\l2\<\k\+n^/3 |ii|<|fc|+„i/3 



Hence, using (|2.2Up .we get (|2.42p . Now let us observe that 
y'(A) - Ay"(0) ,^ 1 £ 



f V\X)-XV"{0) ^ n( \ 

J A^ p{X)dX = -^Q{p) 



/i=0 ^ 



Hence, 



= ^(^"(0))' + iYv^'Hj^-'^)^^,,^^,, (.(^)(A) = iv'ix) - xv'mx- 
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where 

c(2)= I V^'HX){pn,nW - piX))dX. 



Using Corollary [2] from Lemma [3l we get 

v^'\j^^^)n+k,n+k = I ^^^=§- + O m/n) + 0(n-2/3) = Po(0) + O i\k\/n) + 0(n-2/3). 
Therefore 

4'i = ^(^"(0))^ + ^o(O)^ + O {e/n') + c(2). (3.38) 

Now we apply ()3.10p to 

(,) _ V'iX)-XV"{0)-lX^V(^\0) 
^ ' ^" X^{z-X) 

We get 

|fc| 



.(4) 



/ (A, z)p(A)dA ± n-i j("))„±i,n±i + O (n-^^ log^/^ n/|9z| 

= ^^^M + o(z) + O iQgi/2 ^/iQ^i") + o {k/n) 



= ^ (24Po^(0) + 2y'(0)y"(0)) + 0(z) + O (n-^/"^ log^/2 n/|9z|) + O {k/n) , (3.39) 
where the last equality follows from (j3.36p . Now, substituting (j3.37p - (|3.39p in p.SSp . we find 
1 / F^^^foA / k ^ 

5„+fc,„(z) = - ZV"{0) + ^ - ( -^0 (0)^' - -^O(O)^' -Ck + 6n+k,n{^) - Sn+k,n{z) 



(3.40) 

with Cfc defined by (j2.42p . 6n+k,n{z) defined by (j3.34p and 

W(^) = 4°^ + ^'4') + (0(n-^/3) + 0{k'/n')) + 0(^5) + ^^Oin-'/' logV^ „). (3.41) 

Since p.4ip follows from (I3.40p . we are left to estimate Cn \ Cn^ and Sn+k,n{z). 

Taking into account (j3.34p . to estimate Sn+k,n{z) we need to estimate {R'^)n+k,n+k and 
{R'^)n+k,n+k-i- Let us take N' = k + log^ nn^^^, and consider J{N') defined by (j2.3p and 

i?«(z) = {z- J^"") - J{N')y\ 

Then, using the resolvent identity (j2.64p and (j2.62p . we get for any z : 9z > n^^/^ 

\{RX+k,n+k - {R^'^ * R^'\k{z)\ < Ce-^'°s^". (3.42) 
Applying the resolvent identity (12.64p to R^^^ = (z — J'^^^)^^ and R^^\z) defined above, we get 

|(i?«*i?«)fc,fc(z)-(i?W*i?W)fc,fc(z)| 

< 2 ((ii(°n^) * R'^'Hm^k + \Qz\-\R<^^\z) * J\N') * i?(°)(z)),,fc) 

< ^ + ^ + ^ E - J1 < f 1 + 4^ + ^) • (3-43) 

\ksz\ n\ksz\-^ n\^z\'^ ^-^ \<sz\ \ n|yzp n|ksz|'^y 
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Here we have used that for |92:| < 1 



xsz\ 



Substituting (fOHj) in (lOD . we get the first estimate in (ITIHD . 

To estimate c>^ and c^f'* we subtract from (j3.33p the same equation for k := k — 1 and 
multiply the result by n (see the proof of Lemma [T] for the details) . Then we get 

2gn+k,n{^)Rn+k,n+k{z) = J ^^H^^nlkWfdX - (3-44) 

where 

n+k 



<^i+fe,n(^) = - ^Rn+k,j{z)Rj^n+k{z) - Rn+k,j{z)Rj,n+k{z)- 

Using the same trick as in p.42p . we get 

oo k 

Besides, since ii[,'^j(2;) is an even function of (j — k), we observe that 

j=—oo j=—oo 
( R) 

Hence, to estimate ^(-z) it is enough to estimate the difference between r.h.s. of the last 
two formulas. Using (j3.43p for the difference of the first sums, similar bound for the difference 



of the second sums, and the bound = — 4| ^/"^ < C for \z\ < 1, we get 



ie,„WI<;|^(l + + (3.45) 

Now performing transformations ()3.35p for the integral in the r.h.s. of ()3.44p . we can rewrite it 

as 

Rn+k,n+k{z){2gn+k,n{^) " ^"(0)^ " V(^\o)z^ /6) = af^z^ - afl + 6i\„{z) + Oiz"), (3.46) 
where 

4°i = v('\j(^^)Uk,n+k = 2n-2/3Po(0)g2(-^) + |- + 0(^-1) 

«S = v(^\j(^^)n+k,n+k = Po(0) + 0(n-2/3). 

Let us take k > and change the variable z = eC with = k/PQ{0)n in (j3.46p . Then, using 
([Oo]) . we obtain from ([3^ 

C^ + ^ + 2^0(0)^g^(;^)+6-^P0-'(0)C,n(^'C) + o(l) _ i?,(C) 

^n+k,n+k{£i) 



2i(^C^ + + e-^^V'(O) (cfc + 4+fc,n(eC) - <5n+fc,n(eC)) ) 



(3.47) 
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In view of (f05]) 



^~-'l45(^T)|-0, as 



nV3 



cxo, 9^ > d 



with any fixed d (see (|3.45|) ). Besides, — > 0, as x — > oo, because of (|2.32|) . Therefore there 
exists some fixed Iq > 0, such that for k > Iqu^^^ and any ( : > 1/4 



1 

< — < min 

4 |C-i/V2|=l/4 



Then according to the Rouchet theorem i?i(C) has a root inside the circle B of radius 1/4 
centered at Thus, if R2{C) has no roots of the second order inside B, then similarly to the 

proof of Lemma [T] we obtain a contradiction with (j3.17p . Therefore, using the first inequality of 
()2.43p . ()3.41|) and (I3.45P we conclude that there exists an absolute constant Co, such that 

|cW|,|c^(f)p<Con-^/l 

These bounds and (|3.4ip prove the second estimate of ()2.43p . 



^(A)-S(A))V'it(A)V'it+,(A)dA 



4 Appendix 

Proof of Proposition d Using the spectral theorem and Proposition [71 we get 

n+k^n+k+S 

Let us represent v{X) by its Fourier expansion 

oo 

v{\)= v,e^^'\ d 

j=-oo 

Then we have 



< Ce 



-nCe 



vr 



(4.1) 



v{J 



j=-oo 



\j\<cM 



+ 0{M-'+^/^), 



(4.2) 



where c is some absolute constant which we will choose later. The bound for the remainder term 
in the last formula follows from the estimate 



\j\>cM 



< 



E 

j\>cM 



\Vj\ < 



E 

\j\>cM 



|2| -12^ 



\ 1/2 / \ 1/2 

) ( E i^r'') <ii-^')|i2(cM)-^+v2 



■\j\>cM 



Consider now N' = [N + M] + 1 and denote by ^("''^') the matrix whose entries coincide 
with that of ^^"^ with the only exception J'j:^^^,^^_^^,_^-^ = 0. We will use the identity, valid for 
any matrices J'l , ^2 



Ji)e''^'ds 



(4.3) 
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Let us take |A;| < iV and apply ([331) to Jn%^n+k+5 -^l+kj+k+s- T^^^"^ we get 

viJ^'^'')n+k,n+k+5 — v{J^'^''^^)n+k,n+k+S 
t 



f J ( I iid(t-s)j'-"''^'\ n-i!") I isj^"\ 

dS Vj Ue )n+k,n±N'Jn±N',n±N'+l\^ )n±N' +l,n+k+5 

•^0 |i|<Af ± ^ 



Now we use the bound, valid for any Jacobi matrix J with coefficients Jk,k+i = Jk+i,k = G R, 
\0'k\ ^ ^- Then there exist positive constants Co, Ci,C2, depending on A such that the matrix 
elements of e**"'' satisfy the inequalities: 

\{e''^)k,j\ < Coe"^^l''~^'l+^'*. (4.5) 
This bound follows from the representation 

JtJ\ _ -*- 1 itz 



where R={zl - J)"\ and from ([2:62]) . 

Using ([I3D in ([0|l . we get for any c < Co(Cid)"^ 

^(j("^)n+fc,n+fc+5 - v{J^^''''\^u,n+k+6 = 0{{cM)-'+^'^) . 

Similarly (see definitions ([2.2p and ([2.3p ): 

t;(j(°) + J)k,k+s - v{J^''^'''\+k,n+k+5 = 0{{cM)-'+'/^). 



These two bounds give us ([2.4p . 

Proof of Proposition\^ Assume that \xk\ > £ for some k. Without loss of generality we can 
assume that > 0. Then due to ([2.33P 

Xk+i - 2xk + Xk-i > xl. 

Consider first the case when also = Xk+i — Xk > 0. Then by induction for any M — k > i > 
we have > d^^\ Xk+i > xj, and > Hence 

ei > Xfc+Mi >Xk + xlMl/2. 

If < 0, then according to (1233]) we have x_fc_i > x^ and we obtain ([2.34p moving from k 
in the negative direction. 

Similarly, assume that at some point \k\ < M — 2Mi 

> 4max{e^ (2eiAff 2)2/3} 4^2_ ^^ g^ 

Since |4^^l — ^Ai^ because of ([2.33P and the first inequality of (|2.34p . we have for any |i| < zq := 

[2/(3^)] < [2/(3£)] < Ml 

4'i > 4'V2 ^ > \xk+s^o\ > *o4'V2 = 3/i-i4'\ 



29 



where s = signx^. The last inequahty here contradicts to (14. 6p . Hence, ()4.6p is false and we 
obtain the second inequality of (j2.34p . 

To prove (I2.36|) observe that if we consider two (2M + 1) x (2M + 1) Jacobi matrices J^-^^ 
and J^"^) with entries 

j(/)=Z)(/)-jW, jW = Z?W-jWz)g= 4,^.(2 + /,), Dlf^ = 5,,i2 + d% \k\<2M 
then, using the Neumann expansion for their inverse we get 

< (J(^))^] < (J(^))fc] < (2d)-ie-'^l^-^l. (4.7) 

Hence, rewriting ()2.35p as 

( J^-^^x)fc = -rk + Sk,2MX2M+l + h-2MX-2M-l 

we get (lOHl) from (lOl . 

Proof of Proposition [7| It is evident that it is enough to prove (I2.52p for the case when 
2Po(0) = 1 in (jl.lTp and (j2.52p . Hence, below we consider this case. For x > the statement is 
evident. Let f{x) = y^— x/G — q{x) and xq be the first negative root of /. Since it is known (see 
[12] ) that q{x) = Ai{x){l + o(l)) as x ^ +oo, we conclude that ^(0) > Ai{<S) = 0.355028... > 
3^/^/6 (for Ai{0) see pj) Besides, it is known that q{x) > 0, q'{x) < (see [T2]) and so 

V-xo/G = q{xo) > q{0) > 32/3/6 ^ xo < -31/3/2- 
But for any point x < xq < —3^^'^/2 in which q{x) > y^— x/6 

q"{x) < y/-x/6{x - x/3) < -^^-x3/6 < - (4V6(-xo)3) < (7^^)" . 

Therefore /"(x) > for x > xq. Since by definition /(xq) = 0, /'(xq) < (because /(O) < an 
xo is the first root of /) we conclude that /(x) < for any x < xq that contradicts to (jl.lSp . 
Thus we have proved that the left hand side of (j2.52p is always positive. But since it tends to 
infinity as x ^ ±00, we conclude that there exists positive 5, satisfying ()2.52p . 

Proof of Proposition\^ We use the estimate for matrix elements of the resolvent of an arbitrary 
Jacobi matrix J , with entries < A: 

\T^kM\ < ^e-^^l^^ll^-^-l, (4.8) 
\<sz\ 

This estimate is similar to well-known Combes- Thomas estimates for Schrodinger operator (see 
e.g-liOj). 

Let be the Jacobi matrix, whose entries coincide with that for J' with the only 

exceptions ^^^^^^±^+1 = ji±M+i,k±M = 0' and 7^('='*^) = {z - ji^^^^yK Then, by the 
resolvent identity (j2.64p 

T^k,j — "^fe^/^^ = ^ (j^\^^k±MJk±M,k±M+lT^k±M+l,j + '^l^^k±M+iJk±PI+l,k±MT^k±M,j^ ■ 

Since J^^'^^ has a block structure, its resolvent T?.^'^'*^-* also has a block structure and its 
coefficients Ti^^f^^ depend on Jj^j+i with \j — A;| > M. Hence, we can apply (|4.8p to 

'J^^kj^\ '^lt,k±M a-nd T^^k^'i^^M^i- Then we get (j2.62p . Proposition [6] is proved. 
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